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Abstract 

Field-dependent maxima in the relaxation rate of the magnetic molecules 
Mni2-Ac and Feg-tacn have commonly been ascribed to some resonant tun- 
neling phenomenon. We argue instead that the relaxation can be understood 
as purely due to phonons. The rate maxima arise because of a Jahn- Teller- 
like distortion caused by the coupling of lattice phonons to degenerate Zeeman 
levels of the molecule at the top of the barrier. The binding energy of the 
distorted intermediate states lowers the barrier height and increases the re- 
laxation rate. A nonperturbative calculation of this effect is carried out for a 
model system. An approximate result for the field variation near a maximum 
is found to agree reasonably with experiment. 
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In the last five years, a remarkable and novel type of paramagnetic relaxation has been 
observed in certain molecular crystals (H-Q. The best known, a Mn 12 -acetate complex (or 
just Mn 12 for short), has a ground state spin S = 10, and a four-fold symmetry axis z, which 
coincides with the c axis of the tetragonal crystal that it forms. It is found that if the spin 
is first polarized along c, and then allowed to relax by reversing the field to a value H, the 
relaxation rate is a non monotonic function of H, showing maxima at roughly equally spaced 
fields H n = nH*, with H* = 0.44 T. Hysteresis loops measured with a slow field sweep show 
sharp jumps in the magnetization M at the same reverse fields H n . The relaxation rate V 
obeys an Arrhenius law with a barrier of ~ 62 K when H = 0. Similar effects have been 
seen in [(tacn) 6 Fe 8 2 (OH) 1 2] 8 +, though not as sharply. 

In this Letter we contend that at present there is not even a qualitative understanding of 
these rate maxima, and attempt to fill this gap. We begin by reviewing the basic background 
theory, which is due to Villain, Hartmann-Boutron, Sessoli, and Rettori (VHSR) ||. The 
principle terms in the spin Hamiltonian are 

H s = -DS 2 - hS z , (1) 

where h = g/iBH. The next order terms, and S±, are small and have been ignored. The 
spin states \m), where S z \m) = m\m), have energies e m = —Dm 2 — hm. Suppose the initial 
spin state is | — S), and e_s > eg. As explained by VHSR, the spin relaxes via a multistep 
Orbach process. It absorbs phonons and climbs from m = —S to —S + 1, —S + 2, . . ., until 
it reaches the highest level m*, and then descends to levels m* + 1, m* + 2, . . ., S by emitting 
phonons. The occupation probabilities {p m } obey a set of master equations 

' // '"' -ET> m + Efc (2) 



dt 



n 



where 7^ is the rate or probability per unit time for the m — > n transition. To calculate 
{7™}, VHSR use Fermi's golden rule with a standard spin-phonon interaction 

^s P = 9ijkl u~ ll2 e ai q aj (a a + al)(S k Si + SiS k ) , (3) 

a 

where a labels the phonon modes with creation and anihilation operators and a a , q a , e a , 
and oj a are the wavevector, polarization, and frequency, and i, j, etc. are Cartesian indices. 
The tensor is only weakly q dependent. With e mn = e m — e n , we obtain 

rJn _ tr\\/ 1 2 v J e nm n ( e nm)i e m < e n> / A\ 

Tm-^l^mnl * S 3 m , / 



/ 1 ■ 



Here, if is a constant dependent on the sound speed, the density of the solid, and the 
couplings Qijku Vmn denotes a spin matrix element such as (m\SkSi\n) , and n(e) is the Bose 
function n(e) = (e^ e — l) -1 , with (3 = 1/ksT. VHSR solve these rate equations and find 
that for H = 0, 

T = K\V 10 \ 2 D 3 exp(-(3DS 2 )/(l - e^ D ). (5) 

The factor e~^ Ds2 yields the expected activated behavior with an energy barrier DS 2 . The 
spin phonon interaction required to fit the observed rate is somewhat large, but not unrea- 
sonable. 
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Since the relaxation is a sequential process in this picture, T is limited by the weakest 
link in the chain j7| . This leads to a disaster whenever the field is such that the topmost two 
levels, m and m' , say, are degenerate. For then, by Eq. (f|), 7™, and 7™ vanish as /c^Te^/. 
One link in the chain is broken, and so T must vanish completely! 

In fact this conclusion is not correct in detail. To see this, we first note that the S X S Z 
and S y S z terms in Eq. (|3]) connect levels with Am = 1, while S x S y like terms allow Am = 2 
processes. (The Am = terms lead to ignorable level shifts.) Second, there are two types of 
level crossings: m and m + 1 coincide (Am = 1 or "odd" crossings) when h = ft,_( 2TO+ i), and 
m and m + 2 coincide (Am = 2 or "even" crossings) when h = ft,_( 2rn+ 2)> where hj = jD. 
Since the Am = 2 channel is open at a Am = 1 crossing, and vice versa, relaxation is still 
possible. Nevertheless, T should be a minimum whenever h = h n for any n, even or odd, as 
one channel is closed off. This still poses a problem, as these are just the field values where 
the observed T is a maximum. Taking D = 0.57 K from electron spin resonance (ESR) 
experiments ||, we get Hj = 0.42j T, which agrees quite well with the measured jump fields 
in the hysteresis loops. 

Given the above agreement, it is extremely tempting to argue that there is a resonant 
coupling between degenerate energy levels near the top of the barrier || due to some pertur- 
bation, such as hyperfine or dipole-dipole interaction, or terms such as S+ + St in the spin 
Hamiltonian |5,10]. By and large, all these mechanisms are either too weak to explain the 



observed widths of 0.03-0.05 T in the relaxation rate maxima, and/or run into difficulties 
with selection rules. The hyperfine field, e.g., has an estimated rms strength of 0.015 T, 
and must act in second order to connect levels with Am = 2, which implies a width ~ 10 -3 
T. Similarly, the dipolar field is estimated as ~ 0.01 T. Not only is this again a little small, 



but an experiment on a frozen solution of Mn 12 |l 1] with greater intermolecular separation 
than the crystal shows faster relaxation, exactly the opposite of what would be expected. 
The S± terms cannot explain the Am = 1 maxima. The strongest evidence that a perturba- 
tion involving only the spin is not responsible for the rate maxima comes from experiments 



in which the field is applied at an angle 9 to the z axis [12,13 1 . Since the perturbation 



—h x S x oc 9 for small 9, T should rise dramatically as 9 increasess from say 0.5° to 10° [14 



with H z held fixed at H n . However, as can be seen from Fig. 4 of Ref. ]nj, the jump AM 
hardly changes at all for H z = Hi, and increases only by 2 at H z = H 2 as 9 is varied from 
to 10°. Further, the width of the jump at Hi is essentially unchanged even as 9 is increased 
to 44°, whereas a simple anticrossing would yield a width varying as .Dsin#. 

How, then, does the magnetization relax when H = H n 7 In our view, the causative 
agency is still the spin-phonon interaction but it now plays a dual role. Specifically, the 
part of TC sp that couples two levels can not be treated perturbatively when these two levels 
are degenerate. One way to see this is to map the degenerate levels onto the s z = ±1/2 
levels of a pseudospin with s = 1/2, to write H sp and TC P as c a s x x a and w a (p„ + x 2 a )/2, 
where x a and p a are oscillator position and momentum coordinates for mode a. It is then 
obvious that the lowest energy eigenstates are those of s x , not s z , and that in these states 
the oscillators have a non zero mean displacement, (x a ) = —c a s x /u a . The binding energy 
Q due to this displacement, effectively leads to a reduction of the net barrier, and hence to 
an enhanced relaxation rate. 
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FIG. 1. 

Upper four magnetic levels in the degenerate case, without (solid lines) and with (dashed lines) 
lattice distortion. 

In the rest of the Letter we shall sketch a calculation of this effect for a model problem, in 
which the spin-phonon interaction is not the most general one possible, but which would also 
yield a vanishing rate at degeneracy if treated perturbatively. We will show that phonons 
alone give a nonzero rate, which in fact exceeds that which would be guessed from a Arrhenius 
law which the barrier varies smoothly with H and the prefactor is taken as constant. We 
will also show that the rate is a local maximum at degeneracy, with an approximate field 
dependence that describes the data on Mni 2 fairly well. 

Our model focuses on the four energy levels near the top of the barrier (see Fig. 1), which 
we label 1, . . ., 4. The full Hamiltonian 7i is the sum of spin, spin-phonon, and phonon 
terms: 

W. = E^L7)01. (6) 
i=i 

Wp = £^(P& + 4), (8) 

where Uj = and we have set h = 1. We will take e\ « €4, e 2 ~ 63, and make the 

inessential simplification a 34 = an. Also, we choose a spectral density 

J ( u ) =nE^( u ~ w ") = Bu 3 exp(-u/u} c ), (9) 

which describes a coupling to accoustic phonons as per by Eq. (0). The cutoff uj c is like a 
Debye frequency. Defining D = e 2 i, we assume u c ^> D, ksT and also that the maximal 
lattice distortion binding energy Q (realized in the case of exact degeneracy e 2 = £3) obeys 
fiCD. This energy is given by 

n = ^fd^ = c&E-^-- (io) 

71 Jo u 23 Y2w« v ; 

Note that since all oscillators contribute to Q and not just those near uj = 0, the high 
frequency details of J(oo) can significantly affect Q. 

We now wish to calculate 7*, the net rate at which the spin, initially in state |1), relaxes 
to state |4), at a temperature T. For future use and for comparison with Eq. (|j), we record 
here the result of a Fermi golden rule calculation of the 1 — > 2 transition rate: 
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7a ° bs = 2a\ 2 J{D)n{D) » 2Baj 2 D 3 / (e? D - 1). (11) 

A similar golden rule calculation gives 7! = 0, and as discussed earlier, utterly fails to 
explain how the spin relaxes. The correct procedure is to first find the exact eigenstates of 
the full Hamiltonian in the 23 subspace, H23 = -P23^-P23, -P23 = |2)(2| + |3)(3|, and then find 
the absorption and emission rates to and from these intermediate states by the golden rule. 
The absorption rate from state 1 so calculated can be formally written as 

7ab s = 27r^^(l,//|y t 5(e 1 + e M -W 2 3)y|l ) //), (12) 

where V = a±2 Y, a 

c a x a |2)(l|, and fi labels the initial phonon states, with energy 6^, and 
occurrence probability p M cx exp(— /3e M ). The emission rate can be similarly expressed. 

Formula ( |12"D can be greatly simplified when e 2 = €3. Fourier transforming the energy 
^-function, and expressing the sum over /x in terms of the density matrix exp(— /3TC P ), we 
obtain a trace over phonons which involves only Gaussian integrals. A lengthy but simple 
analysis yields 



r d re-^-^ R ^ U(r) + 4,S*(r)l , (13) 
J —00 L J 



where 



Tabs — a l2 



S 1 (t) = - [°°dajJ(u)(f a + f e ), (14) 

7T JO 

S 2 (T) = - fdu;^(l + f a ~ f e ), (15) 
7T JO UJ 



1 /- 00 , J(u) 



1 r°° 
r( t ) = - du 
n Jo 



uU 2 



ujt + i coth ^ - i(f a + f e 



(16) 



with f a = e iulT n(u), and f e = e~ iu)T [l + n(u)]. 

The essential aspects of the integral (IT^) follow from the behavior of a\^R{r) as r — > 00. 
The first term in Eq. ([IB]) grows linearly with r, the second tends to a constant, and the 
last two vanish as r — > 00. We make this explicit by writing 

R( T ) = a^iyir + iW) - iR'(r), (17) 

where Q is the binding energy ([TOD- 



W = ^ rfcjA^coth— ,. 18 
ix Jo LU 2 2 



and R'(t) = f du)J(u)(f a + f e )/ixu) 2 . The term Qr clearly must be combined with Dt in 
the exponential factor in Eq. flT3"|). This point makes it obvious that the relevant transition 
energy when £32 = is not the bare D, but rather the reduced quantity 

D = D-a (19) 

Likewise, the term W reduces the coupling constant a 12 to a 12 = a i2 exp(— W/2). This is 
nothing but the Franck-Condon overlap factor between the unpolarized and polarized lattice 
states accompanying the spin states 1 and 2 (or 3), respectively. 
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The explicit form (H) for J(uo) yields f2 = 2Bal 3 LU 3 /ir, an d W « Bal 3 uj 2 /n = 2Vl/uo c for 
fc^T <C w c . Since £7 < D and D u c by assumption, W <C 1. 

We can further approximate Eq. ( |i~3"D by expanding the remaining factor expfa^-R^ 7 ")] 
in powers of a 2 3 up to 0(a| 3 ), which is justified as long as a 2 3 is not too large, and reflects 
what we mean by the "weak" non-perturbative limit. If this is done, 7 a b s can be written as 
the sum of "direct" and "indirect" parts 



/oo ~ . . 

dre-^S^r) = 2a{ 2 J(D)n(D), (20) 
-00 

dTe~ WT \S^r) + S 1 (r)R'(T)]. (21) 
-00 L - 1 



7in — a l2 a 23 



The last result for 7di r follows by first doing the r integral, and noting that D > 0, so only 
the absorption term f a in Eq. (H) contributes. A similar procedure shows that 7i n can be 



expressed as the sum of three parts, 

7in,i = 4a? 2 £j(^)np), (22) 

7in,2 = -ct>l 2 al 3 D 2 [ [ du) x du)2 Ul f w(a>i)n(a> 2 )6(a>i + uo 2 - D), (23) 

7in3 = -a^a 2 ^ 2 f[du Jl du J2 Ji - UJl l J r ^ 2) n(a;i)[l + n(a; 2 )]5(a; 1 - w 2 - D). (24) 

It can now be shown that 7^,2/7^,1 ~ (D/u c ) 3 for fc^T <C D, and ~ D 2 knT/uj 3 for /c^T ^ 
D, while 7in,3/7in,i ~ D{k B T) 2 /ui 3 for all T. It is assumed that k B T <C u c always. Dropping 
7 in 2 and 7^3 completely, we obtain the compact approximation [compare with Eq. ([TT|)1, 

7abs « 2a 2 2 (l + 2tt/D)j(D)n(D) 

= 2Ba 2 12 D 3 (l + 2n/D)/( e f 35 - 1). (25) 

Note that 7in <C 7 di r because fl « D, which renders the expansion in a 2 3 sensible. The 
condition Q -C D is actually basic to the very formulation of our calculation, for otherwise 
the levels |1) and |4) will also induce significant polarization of the lattice, and treating them 
as bare levels is unjustified. 

The last step is to find the over all rate 7^. Let us do this first for the symmetric 
case, where e 32 = e 4 i = 0. Then, phonon absorption starting from state |1) takes the spin 
to states |r, 1} = 2~ 1 / 2 (|2) ± |3)). The individual rates for these transitions clearly obey 
7i = 7i = 7abs/2. Liekwise the rates for phonon emission from these levels to level 4 are 
equal, and given by 7^ = 7/ = e /3D, y^ s /2 using detailed balance. Following Ref. ||, one 
finds 

7i = 7abs- (26) 

We thus see that even with degenerate top levels, phonons can lead to spin relaxation. 
In fact the rate (EBf) exceeds that given by a classical Arrhenius law in which the field 



dependence is smooth and enters only through the overall energy barrier. In making this 
comparison, let us also incorporate the lower 2S — 3 energy levels into the picture. The net 
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rate is again determined by the slowest step, and at low temperatures, exp(— (3D) <C 1, it 
can be found by multiplying 7^ by the Boltzmann factor for level 1. Ignoring factors of order 
unity, and writing the matrix element as \V\ 2 , we have 

T ~ K\V\ 2 D 2 (D + 2n)e-W+W-n)l } (27) 

where U — e\ — is the full energy barrier. This should be compared with a "naive" rate, 
r° ~ K\V\ 2 D 3 e~P u . It is not hard to see that T > T° provided k B T < min(cj c , D 2 /Vl) [TjJ. 

We can also ask for the rate in the asymmetric case, when €32 7^ 0. We do not now 
have closed form answers, and our results are more approximate. Let €2 = —63 = C/2, and 
€\ = —€4 = (, where ( <C D. Lattice polarization will now change ( to (e~ w , and the spin 
states will be slightly rotated from |r, I). The main effect, however, will be that the binding 
energy itself is reduced: Q — > Q — c( 2 /Q, where c = 0(1). This means that the center of 
gravity of the levels 2 and 3 will not be lowered as much as in the degenerate case. Secondly, 
we can safely approximate the rates 7^, etc. by the form (|25| ) (times 1/2), but with the 
correct energy differences in the Bose factors. The rate 7^ can now be written as a sum of 
rates through the parallel channels |r) and |/). A straightforward analysis shows that the 
terms linear in ( (due to differences in 7]* and r y[, e.g.) cancel. The chief surviving term is 
from the systematic reduction in binding energy, leading to a factor exp(— c(3( 2 /Q) in 7^. 
This shows that the relaxation rate has a local maximum near the coincidence fields ||17|| . 

We conclude by asking how far our model applies to real Mni2- The obvious omissions 
are (a) Am = 2 processes, and (b) tj^ — tkj couplings in Eq. (0). Neither of these is expected 
to change the qualitative picture given here. The key unknowns are the binding energy Q, 
and whether Q < I). A look at Fig. 2 of Ref. |T3j shows that a field variation of the form 
T(H) ~ exp(— c/3( 2 /Q) with ( = Am(h — h n ), describes the rate maxima quite well, and we 
find fl/c = 9, 6, 10, and 4 mK for n = 0, 1, 2, and 3, respectively. Direct estimates based, 
e.g., on Eq. (B.3) of Ref. [i~0||, yield Q values 0.5-20 mK, which are quite consistent given 
the uncertainties in uj c etc. Direct experimental probes of 7i sp are clearly very important to 
testing our theory. Suggestions for how this may be done, and more detailed calculations 
will be published elsewhere. 

I am indebted to B. Barbara, A. J. Leggett, C. Paulsen, M. Randeria, and R. Sessoli for 
many useful discussions. This work was done in part at the Institute for Nuclear Theory, 
University of Washington, Seattle, and I thank the memebrs and staff of the Institute for 
their warm hospitality. This work is supported by the National Science Foundation through 
grant number DMR-9616749. 



7 



REFERENCES 



* Electronic address: agarg@nwu.edu 
[1] M. A. Novak and R. Sessoli, in Quantun Tunneling of Magnetization — QTM'94, edited 

by L. Gunther and B. Barbara (Kluwer, Dordecht, 1995). 
[2] J. Friedman, M. P. Sarachik, J. Tejada, and R. Ziolo, Phys. Rev. Lett. 76, 3830 (1996). 
[3] L. Thomas et al, Nature 383, 145 (1996). 
[4] C. Sangregorio et al, Phys. Rev. Lett. 78, 4645 (1997). 

[5] J. Villain, F. Hartman-Boutron, R. Sessoli, and A. Rettori, Europhys. Lett. 27, 159 
(1994). 

[6] A. L. Barre, D. Gatteschi, and R. Sessoli, Phys. Rev. B 56, 8192 (1997). 

[7] Thus, Eq. (|j) can be understood as the Boltzmann factor for occupation of the level 
m — — 1, times the rate 7° 1 . 

[8] The simultaneous degeneracy of lower pairs of levels when h = nD implied by Eq. ([I]) 
is surely artificial. ESR indicates a correction term B(S^ — S?) for Fe 8 with B = 0.046 
K, and D'S* for Mn 12 with D' = —1.1 x 10 -3 K. These lead to nondegeneracies for the 
lower levels that are comparable or bigger than the perturbations used to couple these 
levels. For Mni2, e.g., we obtain E\ — E-2 = 12D' when h = hi, and E\ — E_3 = A8D' 
when h — h,2- 

[9] A. L. Burin, N. V. Prokof'ev, and P. C. E. Stamp, Phys. Rev. Lett. 76, 3040 (1996); 
A. Fort et al, ibid 80, 612 (1998); D. A. Garanin and E. M. Chudnovsky, Phys. Rev. B 
56, 11,102 (1997). 

[10] F. Hartmann-Boutron, P. Politi, and J. Villain, Int. J. Mod. Phys. B 10, 2577 (1996). 
[11] R. Sessoli, Mol. Cryst. Liq. Cryst. 274, 145 (1995). 
[12] J. R. Friedman et al., J. Appl. Phys. 81, 3978 (1997). 
[13] F. Lionti et al., J. Appl. Phys. 81, 4608 (1997). 

[14] This argument would fail if the relaxation was already saturated at 0.5° (A. J. Leggett, 
comment at seminar). This is unlikely to be so, because the jump at H\ is very small, 
and less than 3% of the total AM for all 9 < 10°. 

[15] The names are based on the prefactors in Ojj+i, and should not be interpreted too 
closely. 

[16] It is quite likely that the "naive" rate T°(H) can not quantitatively describe Mn 12 for 
any value of H at low temperatures. The comparison made here serves only to prove 
the point that lattice distortion effects are large. 

[17] At very low temperatures, the systematic increase in from the decrease in e 2 i may 
nullify this effect, but we do not have any hard calculations to support this point. 



8 



